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C3 ABSTRACT 

We place the most robust constraint to date on the scale of the turnover in the 
cosmological matter power spectrum using data from the WiggleZ Dark Energy 
Survey. We find this feature to lie at a scale of fc =0.0160to;(x)35 [ft/Mpc] (68% 
ly-j confidence) for an effective redshift of z c ff=0.62 and obtain from this the first- 

ever turnover-derived distance and cosmology constraints: a measure of the cosmic 
V,Q distance-redshift relation in units of the horizon scale at the redshift of radiation- 

matter equality (m) of I)v(2efi = 0-62)/rH=18.3+3; 3 and, assuming a prior on 
the number of extra relativistic degrees of freedom iVeff=3, constraints on the 
cosmological matter density parameter J7m^ 2 =0.136^q'q52 and on the redshift of 
£NJ matter-radiation equality z cq =3274^2g . We stress that these results are obtained 

t-H within the theoretical framework of Gaussian primordial fluctuations and linear 

large-scale bias. With this caveat, all results are in excellent agreement with the 
• i-h predictions of standard ACDM models. Our constraints on the logarithmic slope 

of the power spectrum on scales larger than the turnover is bounded in the lower 
limit with values only as low as —1 allowed, with the prediction of P(k)ock from 
standard ACDM models easily accommodated by our results. Lastly, we generate 
forecasts to estimate the achievable precision of future surveys at constraining 
ho, f^M^ 2 , z oq and N e ff. We find that the Baryon Oscillation Spectroscopic Sur- 
vey (BOSS) should substantially improve upon the WiggleZ turnover constraint, 
reaching a precision on ko of ±9% (68% confidence), translating to precisions 
on riMh 2 and z eq of ±10% (assuming a prior JV e g=3) and on iV e ff of 1^% (as- 
suming a prior J7m^ 2 =0.135). This represents sufficient precision to sharpen the 
constraints on N e e from WMAP, particularly in its upper limit. For Euclid, we 
find corresponding attainable precisions on (kg, ^m^ 2 , N e s ) of (3,4, l^D^o- This 
represents a precision approaching our forecasts for the Planck Surveyor. 

Key words: surveys, large-scale structure, cosmological parameters 



C3 



© 2012 RAS 



2 Poole et al. 



1 INTRODUCTION 

Recent decades have witnessed an incredible refinement 
of our cosmological model with significant advancements 
in the volume and redshift coverage of galaxy surveys, 
and in the methods of their analysis, playing a pivotal 
role. Such advances have opened up entirely new avenues 
of investigation, notably the use of features in the distri- 
bution of galaxies as standard rulers for direct geometric 
mapping of the Universe's expansion history. 

The most investigated of these features are those in- 
duced by "Baryon Acoustic Oscillations" (BAOs) in the 
early-Universe's photon-baryon fluid. The precise mea- 
surement across cosmic time of the BAO scale in the 
distribution of galaxies - a weak harmonic ripple with a 
comoving fundamental scale of ~105 [Mpc/fr] - has been 
a major driver of recent and future galaxy redshift sur- 
veys and has received a great deal of attention in recent 



literature (e.g. |Cole et al. |2005 


Eisenstein et al. |2005 


Blake et al.|2011b |Beutler et al. 


2011||Blake et al.|2011c 


Anderson et al.|2012 


1. 



Far less studied, the "turnover" in the galaxy power 
spectrum also represents a standard-ruler in the observed 
distribution of galaxies. This feature, which is predicted 
to manifest on comoving scales of ~400 [Mpc//i] (or cor- 
responding wavenumber fco~0.016 [ft/Mpc]), was estab- 
lished at the epoch when the dynamics of the Universe 
transitioned from being dominated by relativistic mate- 
rial (photons and neutrinos) to being dominated by mat- 
ter (dark and baryonic). This transition occurred because 
the mass-energy density of relativistic and non-relativistic 
materials decline differently with a, the expansion factor 
of the Universe. The particle density of both decline as 
p oc a~ 3 but the per-particle energy of relativistic mate- 
rial also decreases with a due to redshift effects. 

Previous to this transition from radiation domina- 
tion to matter domination (at redshift z eq ), oscillations 
smaller than the horizon in the strongly coupled matter- 
radiation field were suppressed by the effects of radiation 
pressure while causally disconnected fluctuations larger 
than the horizon collapsed unhindered (see |Eisenstein fc] 
Hu 1998). As a result, the scale-free primordial matter 



power spectrum believed to have emerged from inflation, 
P{k)ock" s with n s ~l, became distorted such that P{k) 
became a decreasing function of k at small scales with 
a limiting behaviour P{k)cck~ 3 . At the turnover scale, a 
peak in P(k) arose separating scales where P(k) increased 
with k from those which decreased with k. 

While the size of the horizon grew during the epoch 
of radiation domination, the scale of the turnover shifted 
ever larger to smaller values of fco until the epoch of mat- 
ter domination commenced at redshift z cq and the sup- 
pression of small scale fluctuations ceased. From this point 
forward, all scales grew independently and by the same 
fractional amount while in the linear regime, and the co- 
moving scale of the turnover became fixed. The exact size 
of this scale and the time of matter-radiation equality 
was influenced by the matter and radiation mass-energy 
densities; Q.yih 2 and fi r /i 2 , the latter being set by the ef- 
fective number of extra relativistic degrees-of-freendom 
(iV c ff) and by the photon mass-energy density (Q-yh 2 ) 
which is well determined from observations of the tem- 
perature of the Cosmic Microwave Background (CMB). 

Analyses of the CMB have placed the most pow- 
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erful constraints to date on parameters describing the 
epoch of matter-radiation equality by measuring z eq , one 
of the "fundamental observables" of the CMB ( |Komats"u| 
et al.|[2009 |. This measurement is made possible because 



of early integrated Sachs-Wolfe effects which manifest in 
the ratio of power between the 1 st and 3 rd peaks in the 
angular power s pectrum. The latest W MAP results find 
z cq = 3145Ii3 9 ( Komatsu et al.|2011 i. 



Studies of the matter power spectrum on scales of 
the turnover can provide an additional window into phys- 
ical processes operating at (or even prior to) the epoch of 
matter-radiation equality; a time before even the epoch of 
recombination observed by CMB experiments. While the 
precision of the turnover scale as a standard ruler is re- 
duced by the fact that it less sharply defined (fractionally 
speaking) than the BAO scale, its study benefits from ly- 
ing firmly in the linear regime at all redshifts, easing com- 
plications which arise from non-linear structure formation 
and redshift-space distortions. As a result, structure on 
scales of the turnover are sensitive to non- Gaussian pro- 
cesses during inflation permitting interesting new studies 



of early-Universe physics (Feldman et al. 1994 Durrer 



et al. 120031). For instance, scale-dependent bias effects are 



expected for galaxy samples with biases deviating from 
unity 



(e.g. 



Dalai et al. 2008 1, providing a rare opportunity 
to study inflation. Even in the case of Gaussian fluctua- 
tions, informative scale-dependent effects may be present 
on scales beyond the turnover ( |Yoo|20fo| . 

The sensitivity of turnover scales to N B s is also of 
interest given the conflict of several recent CMB stud- 
ies with conventional theoretical expectations. For the 
standard ACDM model, the expected value is jV e ff~3.046 
( Mangano et al.|2005 i.e. slightly larger than 3 owing to 
the fact that neutrino decoupling was not instantaneous) . 
However, recent high spatial resolution CMB measure- 
ments from the ACT and SPT observatories ( |Dunkley| 



et al. 2011 Keisler et al. 20111 suggest a significantly 
higher result of iVeff~4 (Hou et al. 2011 Smith et al.| 



2012 Calabrese et al.|2012 1 leading several researchers to 



explore a class of "dark radiation" models {e.g. Archidia- 
cono et al.|2011 k 



The primary challenge to studies of the turnover, and 
the reason for its scant study to date, is the large volumes 
one must uniformly probe to detect it with any precision. 
While the fluctuations involved may lie on ~400 [Mpc//i] 
scales, many modes on these scales must be enclosed by a 
survey's volume to permit a statistically significant anal- 
ysis. Furthermore, subtle errors in the calibration of a 
survey's selection function across its volume can easily 
lead to significant systematic power spectrum distortions 
on a survey's largest scales. This can be particularly se- 
rious for angular {i.e. 2D rather than 3D) measurements 
in imaging surveys which are much more sensitive to an- 
gular systematics in photometry and/or target selection 
(e.ff. |Ross et al.|20U[[20T2l ). 

Despite these challenges, several attempts have been 
made to measure the turnover or to study the matter 
power-spectrum on its scales. The earliest attempts were 



performed using the APM galaxy survey (Baugh & Ef- 



stathiou 1993 19941 or through studies of the distribu- 



tion of optically-selected galaxy clusters ( Peacock & West 
[19921 |ScarameHa||1993||Einasto et al.||1993| |Tadros et al. 
1998] [Einasto et al.||1999| ). These studies have tended to 
observe turnovers at scales of ~100 — 200 [Mpc//i], which 
is quite discrepant with a wide range of modern cosmolog- 
ical probes and analyses. In all cases, systematic effects 
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induced by survey selection were found to dominate on 
turnover scales or rigorous studies of such issues were not 
performed. More recently, observations of quasars ( [Out- 
ran! et al.|2003 1 and luminous red galaxies (LRGs) in the 



Sloan Digital Sky Surv ey (SPSS; [York et al. 12000 Pad 
manabhan et al.]|2007| ) have managed to measure power 
on scales reaching to those of the turnover but relatively 
little analysis using this information has been performed. 

In this paper we present the most robust measure- 
ment of the turnover scale to date using data from the 



WiggleZ Dark Energy Survey (Drinkwater et al. 20101 



WiggleZ is a large-scale galaxy redshift survey conducted 
with the AAOmega multi-object spectrograph on the 
Anglo- Australian Telescope at Siding Spring Observatory. 
WiggleZ was designed to study the effect of dark energy 
on the Universe's expansion history and on the growth of 
cosmological structures across an unprecedented period of 
cosmic history. The primary science drivers for the sur- 
vey have been the measurement of the BAO scale (e.g. 
Blake et al.||2011b l, of the growth of cosmological struc- 



ture (e.g. Blake et al. 2011a I and measurement of the 
neutrino mass (e.g. |Riemer-S0rensen et al.]|2012[ ) using 
the clustering pattern of UV-selected galaxies, as well as 
studies of the Universe's most actively star forming galax- 
ies out to a redshift of 1.2 
et al.|2012 1 



(e.g. 



Wisnioski ct al. 2011 Li 



In this work we treat the dominant systematic uncer- 
tainty of our study coming from convolution effects of the 
survey selection function and the subsequent extraction 
of the turnover scale in a statistically rigorous way. We 
also perform forecasts for the ongoing and future surveys 
of the Baryon Oscillation Spectroscopic Survey (BOSS; 
Eisenstein et al.|2011| and Euclid ( Laureijs et al.poTT l. 



Our analysis is presented purely within the theo- 
retical framework of Gaussian primordial fluctuations. 
Models which permit non-Gaussian primordial fluctua- 
tions provide for a wide range of large-scale clustering 
behaviours and we will leave their exploration for future 
work. 

In Section [2] we present our method of construct- 
ing a galaxy power spectrum from the WiggleZ dataset 
which is optimised for studies of the turnover. In Section 
[3]we extract the turnover scale from the WiggleZ observa- 
tions and present the resulting distance measurement to 
z=0.62 and constraints on SImIi 2 and z cq . In Section[4]we 
perform forecasts for future surveys (highlighting results 
for BOSS and Euclid), estimating the volume dependance 
of turnover scale measurement precision and of resulting 
cosmology constraints. We also compute the constraints 
in the iVeff— Om/i 2 plane from the CMB observations of 
WMAP and forecast the constraints achievable by the 
Planck Surveyor. Lastly, we present some discussion and 
our conclusions in Section [5] 

Our choice of fiducial cosmology throughout will be 
a standard ACDM model with rj M =0.27, f2 A =0.73 and 
h=0.7, unless otherwise stated. 



2 OBSERVATIONS 

In this section we present our construction of the WiggleZ 
power spectrum for use in our turnover analysis. This will 
include our methods of deconvolving the effects of the sur- 
vey selection function, of coadding the power spectra ob- 
served in seven independent WiggleZ survey regions and 
our treatment of the survey's radial selection function; the 



dominant source of systematic uncertainty in this current 
analysis. 



2.1 The WiggleZ Power Spectrum 

To construct a power spectrum for this study, we followed 

see- 



the general approach described by Blake et al. ( 2010 



Section 3.1). To summarise, we used the optimal weight- 
ing scheme of Feldman et al. ( 1994 FKP), converting red- 
shifts to distances using our fiducial ACDM cosmological 
model. We individually enclosed the survey cones of each 
region within cuboids of sides (L x ,L y ,L z ) and map each 
region's observed galaxy distribution onto a grid with di- 
mensions (n x ,n y ,n z ) using the "nearest grid point" assign- 
ment scheme. Denoting the resulting distribution as n(x), 
we then applied a Fast Fourier Transform to the resulting 
grid, weighted in accordance with the method described 
in FKP (as given by Eqn. 9 of Blake et al. 2010, using a 
weighting factor more appropriate for our turnover analy- 
sis of P = 20000 [(Mpc/h) 3 ] ), to produce n(k). The cor- 
responding values for each WiggleZ region are presented 
in Table [T] In later sections we will require the covariance 
between the bins in our WiggleZ P(k). This was also de- 
termined following the procedure of FKP as presented in 
Section 3.1 (Eqn. 20) of |Blake et aL] ( |2010[ ). We note that 
increasing the characteristic amplitude used in the FKP 
weighting from Po = 2500 [(Mpc/ft) 3 ] used previously for 
WiggleZ BAO studies could amplify systematics originat- 
ing from low-density regions, but in practice we find the 
consequences are not significant for this dataset. 

Both the FKP weighting and the conversion of n(k) 
to our final estimate of the power spectrum, P(k), depend 
critically on the survey selection function, W(x). This 
function expresses the expected mean density of galaxies 
with spectroscopic redshifts at position x, given the an- 
gular and luminosity survey selection criteria. Again, we 



follow the procedure of Blake et al. ( 2010 see Section 2) to 



determine the angular and radial contributions to W(x). 
While angular effects are estimated in precisely the same 
manner as described in Blake et al. (20101 - accounting 



for effects such as the coverage mask, spatial variations in 
dust extinction and radial completeness of each individual 
2df pointing contributing to the survey - our estimation 
of radial contributions to the survey selection function 
required some modification which is described below in 
Section [2~L3l 

While aliasing from the assignment scheme only af- 
fects the power spectrum on scales far smaller than the 
turnover, we nevertheless correct for this effect using the 
method described by Jing ( 2005 ) . Furthermore, since the 
survey target density is obtained from the survey itself, we 
correct for a potential large-scale P(k) bias in a manner 
analogous to the integral constraint of correlation func- 
tions (see Eqn. 25 of Peacock & Nicholson 19911 for ex- 
ample) by applying a boost to the measured P(k) near 
k=Q, generated using the Fourier transform of the win- 
dow function. 

Lastly, to ensure a power spectrum optimised for 
studying scales around the turnover, we have made the 
following three adjustments to the procedure described 



in Blake et al. (20101 



2.1.1 Redshift Range 

The redshift range used in our analysis has two compet- 
ing effects on our analysis. By increasing the range used 
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Table 1. A summary of parameters relevant to our calculation of the WiggleZ galaxy power spectrum for each observed region. 
The values of (L x ,L y ,L z ) specify the dimensions of the cuboid used to enclose each region and (n x ,n y ,n z ) specify the corresponding 
grid dimensions of the cuboid. Volume and iVg a i are the volume and number of galaxies contributing to our measurement in each 
region (spanning the redshift range [z m ini 2 max ] = [0.4, 0.8]) respectively while n is the resulting average number density of galaxies 
in each region. 



Region 


L x [Mpc/fe] 


L y [Mpc/h] 


L z [Mpc/h] 


n x 


fly 




Volume [(Gpc//i) 3 ] 


iV ga i 


n [(Mpc/h)" 3 ] 


9-hour 


899.4 


520.7 


315.5 


256 


128 


64 


0.148 


18978 


1.28 x 10~ 4 


11-hour 


899.5 


520.7 


318.9 


256 


128 


64 


0.149 


20170 


1.35 x 10" 4 


15-hour 


907.7 


694.7 


353.4 


256 


128 


64 


0.223 


30015 


1.35 x 10~ 4 


22-hour 


894.4 


335.1 


338.6 


256 


64 


64 


0.101 


16146 


1.59 x 10~ 4 


1-hour 


891.8 


300.7 


252.4 


256 


64 


64 


0.0677 


8304 


1.23 x 10" 4 


3-hour 


893.1 


313.3 


305.6 


256 


64 


64 


0.0855 


10241 


1.20 x 10~ 4 


0-hour 


891.6 


241.1 


297.8 


256 


64 


64 


0.0640 


7409 


1.16 x 10~ 4 



we increase the effective volume of our measurement, in- 
creasing the number of modes probed by our study at 
scales of the turnover and increasing the signal of our vol- 
ume limited measurement. Unfortunately, doing so also 
leads to a less plane-parallel survey geometry, increasing 
the covariance induced by the survey selection function 
on our largest-scale modes. 

We have performed our analysis on several redshift 
ranges [z min , z max ], varying z min from 0.3 to 0.4 and z max 
from 0.8 to 0.9. While the choice over these ranges does 
not have a strong effect on our results, we have found 
the range [0.4, 0.8] to yield slightly optimal results since it 
maintains minimal covariance between neighbouring P(k) 
bins while providing a near-maximally precise measure- 
ment. Hence, we will use this range throughout in our 
analysis. 



2.1.2 Choice of Fourier Binning 

We perform all our analysis on power spectra binned by 
wave- number (fc). Our choice of binning also introduces 
two competing effects. By increasing the size of our power 
spectra bins, we increase the signal and reduce the covari- 
ance of neighbouring bins, increasing the precision of our 
power spectra. Unfortunately, we are trying to fit to a 
feature at small values of k and the use of bins that are 
too large will prevent us from resolving it. We have experi- 
mented with a variety of bin sizes ranging from 5k = 0.005 
[ft/Mpc] to 8k — 0.02 [ft/Mpc] and have found the use of 
Sk = 0.005 [/i/Mpc] to be optimal. We deemed further 
reduction of this binning unwarranted given the low sig- 
nal at this point for modes larger than the turnover in 
several regions. In all cases, the medians of the |fc| val- 
ues contributing to each bin are used to represent their 
positions. 



2.1.3 Radial Selection Function 

At the mean of our chosen redshift range, the maximum 
scale probed by WiggleZ in directions transverse to the 
line of sight is ~ 500 [Mpc//i] while the comoving distance 
along the line of sight from z m i n =0.4 to z max =0.8 is ~900 
[Mpc//i]. As a result, most of the information in WiggleZ 
on scales of the turnover is contained within radial modes. 
Consequently, our estimate of the WiggleZ radial selection 
function is a significant source of systematic uncertainty 
in our present analysis. 

We have examined several approaches to determin- 
ing the WiggleZ radial selection function based on the 



Summarising briefly: we fit a smooth analytic function to 
the observed redshift distribution N(z) for each observ- 
ing priority band (WiggleZ observations were prioritised 
by magnitude, with apparently-faint galaxies given high- 
est priority and each band representing an equal interval 
of the range 20.0 ^ r ^ 22.5). In this procedure we must 
choose a functional form with which to parameterise the 
observed WiggleZ N(z). In past efforts, this procedure 
has been performed on each of the seven WiggleZ sur- 
vey regions independently using a polynomial of order dy- 
namically chosen to be that above which the reduced-;^ 2 
statistic does not decrease. 

We have reexamined this procedure using several ap- 
proaches and present three here: the standard polynomial 
fit applied to each region independently (i.e. the standard 
approach of past work; denoted the "polynomial N(z) n 
method), a cubic spline similarly fit to each region sepa- 
rately (denoted the "spline N(z) v method) and a polyno- 
mial fit to the total sum of N(z) for all Northern Galactic 
Pole (NGP) fields jointly and all Southern Galactic Pole 
(SGP) fields jointly (denoted the "joint polynomial N(z)" 
method; note, due to differing optical selection, the Wig- 
gleZ NGP and SGP fields have differing radial selection 
functions). The motivation of this last method is to re- 
duce the effects of cosmic variance by combining regions 
but carries the risk of being susceptible to region-to-region 
systematics. The influence of this choice on our resulting 
WiggleZ P(k) is shown in Figure [I] where we show the 
deconvolved WiggleZ power spectrum for the whole sur- 



vey sample (see Section 2.2 below for details). As noted in 



Blake et aL| ( |2010[ ), only results for k < 0.03 [ft/Mpc] are 
significantly affected by the details of our radial selection 
estimation, leaving the scales of the BAO unaffected. 

We have carried our analysis through with each of 
these methods and find relatively little effect on our 
turnover constraints and no qualitative change to the con- 
clusions of this study. The challenge in the procedure of 
determining N(z) for this purpose is to choose a func- 
tional form nuanced enough to capture the selection in- 
duced by the survey strategy and telescope operations but 
smooth enough not to fit to (and hence, remove) the real 
structure we seek to measure. Hence, for this reason we 
favour the joint polynomial N(z) and for all subsequent 
analysis we quote results derived using this approach. 

Lastly, our method includes a full fc-dependent cor- 
rection for misidentified redshifts (referred to as redshift 
"blunders" ) using the survey simulations described in sec- 



tion 3.2 of Blake et al. ( 2010 1. This includes blunders both 



approach presented in Blake et al. (2010 see Section 2.5) 



from constant-multiple shifts due to emission-line mis- 
identifications and a continuous range of mis-identified 
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— Joint Polynomial N(z) 

— Fiducial AC DM P(k) 




k[h/Mpc] 



Figure 1. The deconvolved co-added WiggleZ P(k) using 3 
different methods for estimating the WiggleZ radial selection 
function (a cubic spline or a polynomial fit to each region inde- 
pendently and our preferred method: a polynomial fit jointly 
to the NGP and SGP fields) . Dark shaded regions depict 68% 
confidence intervals and light shaded regions depict 95% con- 
fidence intervals. The solid black line depicts the maximum 
likelihood P(k). The red line depicts a ACDM model incorpo- 
rating a simple model for redshift-space distortions. Residuals 
(in the bottom panel) depict results after subtraction of the 
joint spline maximum likelihood P(k), illustrating the magni- 
tude of systematic effects associated with uncertainties in the 
WiggleZ radial selection function. 



of the starting point of our calculation. We use 10 5 iter- 
ations for each burn-in phase and 2xl0 6 propositions for 
the final integration used to determine the posterior dis- 
tribution of our 50 observed, deconvolved power spectrum 
bins. All chains have been inspected to ensure that they 
are well mixed. To increase the efficiency of our calcula- 
tion, we use the matrix multiplication approach presented 
Eqn. 17) to convolve our P(k) propo- 



Blakeet al. (2010 



sition sets with each region's selection function. We have 
verified that this approach remains equivalent to a full 3D 
convolution on all scales utilised for this current analysis. 

The results of this calculation are presented in Fig- 
ure [l] In this fi gure we see that at scales fc<0.03 [Mpc//i], 
the deconvolved coadded WiggleZ power spectrum places 
useful constraints on the maximum power permitted on 
scales at-and-beyond the turnover, but places little-or-no 
constraint on the minimum permitted power. For com- 
parison, we also show (in red) a ACDM power spec- 
trum in our fiducial cosmology with a simple model for 
redshift-space effects. This is obtained by taking a bi- 
ased non-linear ( |Smith et al.|[2003| i.e. halofit) power 



spectrum from the Boltzmann code CAMB (Lewis et al. 



2000) and applying the redshift-space distortion model 
of Kaiser ( 1987 1 with a Lorentzian damping term (see 



Equation 9 of Blake et al.|20Tla l. The parameter values 
of the galaxy bias (b 2 = 1.18), redshift-space distortion 
parameter (/3=//&=0.69) and the variable damping term 
(cr v =300 km/s) needed for this model were obtained from 
fits to the 2D power spectrum obtained from the WiggleZ 
NGP fields. Figure [TJ illustrates that this model provides 
an excellent fit to the WiggleZ power spectrum over the 
full range of scales presented. 



sky lines into a lognormal simulation including the full 
selection function. The redshift blunder model is based 
on thousands of repeat redshifts obtained during the sur- 
vey. We find this correction to be small compared to the 
P(k) uncertainty for the scales relevant to our turnover 
analysis. 



3 ANALYSIS 

In this section, we will present our extraction of turnover 
information from the WiggleZ power spectrum presented 
in Section [2] as well as the distance and cosmology con- 
straints obtainable from this measurement. 



2.2 Deconvolution of the Survey Selection 
Function 

The WiggleZ survey was conducted within seven separate 
survey regions, each with different selection functions due 
to varying optical selection, region boundary geometries, 
observing conditions and dust attenuations. While much 
of our analysis in later sections will involve direct fits ap- 
plied jointly to the data of these regions, we seek here to 
generate a single deconvolved power spectrum of the full 
survey sample, with the effects of differing selection func- 
tions removed, in order to present the data contained in 
the survey on its largest scales. 

To achieve this we use a Monte Carlo Markov Chain 
(MCMC) approach. Using the Metropolis-Hastings algo- 
rithm, we generate random sets of propositions for the de- 
convolved power spectra we seek, convolve each set seven 
times using the window function of each region and com- 
pute the joint likelihood of these convolved proposition 
sets against the observed power spectrum of each region. 
We compute this likelihood using the full information of 
each region's covariance matrix. We optimise this calcu- 
lation by drawing propositions from a rotated covariance 
matrix constructed from a short secondary burn-in period 
following a primary burn-in designed to erase the memory 



3.1 Measurement of the WiggleZ Turnover 
Scale 

To extract information about the turnover from the ob- 
served WiggleZ power spectrum, we have followed the 



method of Blake & Bridle ( 2005 1. This approach is model 



independent in the sense that it does not take cosmo- 
logical parameters or model power spectra from Boltz- 
mann codes as inputs. Specifically, we fit (using the same 
MCMC machinery and chain lengths descibed in Section 
2.2| the following model, convolved with the WiggleZ se- 
lection functions and jointly fit to all 7 WiggleZ regions: 



log 10 P(k) = 
where x = 



P (1 - ax 2 ) if k < k Q , 



10 J o . 
log 10 P (1 - fa 

In k — In fco 
In fco 



) if k fco 



(1) 



(2) 



We marginalise over Pq, since its interpretation is compli- 
cated by degenerate normalising parameters such as as 
and the bias of WiggleZ galaxies, and over /3 since its 
interpretation is complicated by all the cosmological pa- 
rameters, redshift-space and scale-dependent bias effects 
which influence the overall shape of P(k) on scales smaller 
than the turnover. Hence, we focus here on just two pa- 
rameters: fco and a, which carry information about the 
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Figure 2. Results of fitting our logarithmic parabola turnover 
model (see Eqn. |4j to the observed WiggleZ P(k) as a func- 
tion of the maximum range in k used for the fit (fc max ). Top 
panel: dark shaded regions depict 68% confidence intervals 
and light shaded regions depict 95% confidence intervals. The 
red dashed line depicts the theoretical expectation for a stan- 
dard ACDM model. Bottom panel: the black line depicts the 
reduced-x 2 statistic and the dotted blue line depicts the num- 
ber of degrees of freedom (rirj F). We see from this that our 
measured turnover scale is stable and yields an acceptable fit 
for fc max >0.1. 



Figure 3. The 2D marginalised posterior distribution func- 
tion in the a-ko plane from fitting our logarithmic parabola 
turnover model (see Eqn. |4j to the observed WiggleZ P(k). 
Red contours depict the 68% and 95% confidence regions. The 
white contours depict the loci of several values of n for the 
model P(k^ko) oc k n . From this plot we see that the Wig- 
gleZ P(k) has an inflection at fco~0.016 [/i/Mpc] with n>— 1, 
distinctly different from the asymptotic value for P(fe3>fco) of 
n=— 3. Our fit easily accommodates the theoretical expectation 
of fco=0.016 [/i/Mpc] and n=l for a standard ACDM model. 



position of the peak in P(k) and whether we have de- 
tected a drop in power at scales larger than this peak 
(this is the case if a > 0). All cosmology constraints de- 
rived from our measurement of the turnover will be based 
purely on our measurement of ko . For our fiducial ACDM 
cosmology, the expected value is fco=0.016 [/i/Mpc]. 

Since the real power spectrum has a changing loga- 
rithmic slope with k at scales smaller than the turnover, 
we need to carefully consider the maximum of the range 
of scales over which we perform our fit (denoted fc max ) to 
avoid systematic biases. We have performed our fit using 
a range of fc max and present the results in Figure [2] Here 
we see that the results of our fit for ko are both stable 
and acceptable (as measured by the reduced-x 2 statistic) 
for fc max J>0.1. However, non-linear effects which alter the 
logarithmic slope of the power spectrum are expected to 
become significant for fc>0.2 [/i/Mpc] and so we limit our 
fit with fc max =0.2 [/i/Mpc] for the remainder of our anal- 
ysis. 

The final marginalised results for ko from our fit are 
fco=0.0160±g;ggf| (68% confidence); ig;gj$j (95% confi- 
dence). Our results are not well constrained for a where 
we measure lower limits of a> — 1.32 (95%) with a hard 
prior a<10. In Figure [3] we present the 2D-marginalized 
posterior distribution function of our fit to the WiggleZ 
power spectrum in the a— ko plane. We can see from this 
figure that our analysis constrains ko to lie near the theo- 
retical value of fco~0.016 [/i/Mpc]. Furthermore, a is con- 
strained to a>— 1 indicating that P(k) is either an in- 
creasing function of k or nearly constant at scales larger 
than the turnover. The long tail to positive values of a 
reflects the lack of constraint WiggleZ places on the min- 
imum power allowed on the largest scales of the survey. 

Negative values of a are marginally allowed by our 
fit, permitting a power spectrum which does not formally 
exhibit a turnover. To aid our interpretation of this result 
we have overlayed onto Figure [3] (with white contours) 



a model which presents the logarithmic slope (denoted 
n) for P(k<ko)ock n implied by Equation [4] between ko 
and our largest-scale bin at fc=0.005 [fe/Mpc]. This model 
incorporates the strong and tight degeneracy 



log 10 P = -0.67 log, n fc + 3.12 



(3) 



present in our fit. From this we can see that WiggleZ 
constrains the power at scales larger than the turnover 
to n> — 1 at nearly 95% confidence. The large-fc asymp- 
totic value of n=— 3 is completely ruled out, meaning that 
we have certainly detected an inflection in the logarith- 
mic slope of P(k). With the theoretically expected value 
of n—1 easily accommodated by our fit, and a preferred 
scale for the P(k) inflection of fco~0.016 [/i/Mpc], we thus 
find that our turnover fit is in good agreement with the 
theoretical expectations of a standard ACDM model. 



3.2 The WiggleZ Turnover Distance 
Measurement 

The scale of the turnover roughly corresponds to the scale 
of the horizon at the epoch of matter-radiation equality 
{kn) given by (see Prada et al.|2011 1 

(4) 



where rjj =c 



-2v^)r-H 1 

(1+^r 1 



da 



jo a»H(a) 
with H(a) being the Hubble expansion rate given by 



H 2 (a) = Hq [n r a- 4 +!T2Ma" 3 +fi 



A 



(5) 



(6) 



and Ho being the Hubble constant. This yields a hori- 
zon scale of rH=H7.9 [Mpc] corresponding to /ch=0.014 
[/i/Mpc] for the WMAP result of z cq = 3145 (|Komatsu 



|TaI][20TTl in our fiducial ACDM cosmology. However 
the actual precise position of the peak of P(k) can differ 
significantly from this value. Indeed, using the Boltzmann 
code CAMB ( |Lewis et al.|2000 l we find that the turnover 
scale in our fiducial cosmology is fco,fid=0.0159 [/i/Mpc]. 
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As discussed in Section |2.1.3[ the largest scales of 
the power spectrum utilised in this work are dominated 
by radial modes. However, for most of the scales involved 
in our turnover fit, this is not the case. For this reason, 
we choose to convert our measurement of the turnover 
scale into a distance measurement following the approach 
introduced by Eisenstein et al.| ( |2005[ > for BAO studies. 
Through this approach, our turnover measurement con- 
strains the dilation measure Dv defined as 



Dv(z) = 



-,1/3 



(7) 



where c is the speed of light, Da(z) is the angular diam- 
eter distance at redshift z and H(z) is given by Equation 
[fi] with a — (1 + z) -1 . We will obtain this quantity us- 
ing the assumption that distances scale proportionally by 
the same "stretch factor" (a) under small perturbations 
from a fiducial cosmology. Given the turnover scale in our 
fiducial cosmology, we obtain this stretch factor from 



<S=fco,fid/fco 



(8) 



We will express this dilation measure in a dimension- 
less form using units of the z=3145 horizon size (j"h) as 
dt=Dy /Vh- We choose this scale instead of the turnover 
scale from CAMB to render the result less model- 
dependant. 

For this method we require an effective redshift for 
our measurement, which we compute by determining the 
effective redshift at fc=0.015 [/i/Mpc] for the power spec- 
trum used in this analysis. This is achieved using Equa- 



(2011b) from which we obtain 



tion 13 of Blake et al. 
2cft=0.62; quite similar to the effective redshift of other 
WiggleZ studies. 

The dilation measure at z e g in our fiducial cosmol- 
ogy is dt,fid(z e //=0.62) = 18.41. For our measurement of 
fco=0.0160tgS (68%); la°075 ( 95 %) we obtain a stretch 
factor 5=0.99±g;^ (68%); tlf 7 (95%). Scaling with re- 
spect to our fiducial cosmology using this stretch fac- 
tor, we thus obtain a model independent distance of 
dt(^eff=0.62)=18.3l|;| (68%); ±g 9 7 7 (95%). Assuming the 
fiducial turnover scale given above, this corresponds to 
Dv(zeff=0.62)=2156tS$ (68%); t?g 4 (95%) [Mpc]. 



3.3 Cosmology Constraints from the WiggleZ 
Turnover Measurement 

To calibrate the scale of the turnover as a function of 
cosmology we ran CAMB to produce a series of matter 
power spectra, fixing all parameters to our fiducial cos- 
mology but allowing Qy[h 2 to vary for several choices of 
iVeff. The results are presented in Figure [5] where we illus- 
trate our method of converting measurements of fco into 
cosmological constraints. 

The position of the turnover depends on SIm^ 2 and 
iVcfi in the following way (see Komatsu et al.|2 009). First, 
through the dependance of the total energy density of 
relativistic material on N e g given by 



where fi 



fir =fi 7 [1 + 0.2271AU] 

8ttG 4<7 B Tcmb 



M 2 



(9) 
(10) 



with fi 7 being the energy density of photon^] G New- 
ton's gravitational constant; ob the Stefan-Boltzmann 
constant; and Tcmb the temperature of the CMB, and 




Figure 4. Bottom panel: The ID marginalised posterior dis- 
tribution function (PDF) for fco from fitting our logarithmic 
parabola turnover model (see Eqn. [4| to the observed Wig- 
gleZ P(k). Dark shaded regions depict 68% confidence inter- 
vals and light shaded regions depict 95% confidence intervals. 
Top panel: an illustration of how our constraint on fco maps to 
a constraint on f2jy[/fi r =l + z cq . Lines depict the dependance 
of 1 + z cq on fco, calibrated using CAMB, for several values of 
N e g (increasing from N B ff =0 in magenta to N c g =10 in yellow, 
in steps of 0.67; il^ih 2 increases from left to right). We high- 
light the relation for N e g =3 in red as well as the constraint on 
1 + z cq determined from our measurement of fco for this case. 



second, through the dependance of the redshift of matter- 
radiation equality on fiM and fi r given by 



1 + Z e „ = 



fiM 

~fi7 



fivr^i 



1 



fi 7 /i 2 1 + 0.2271A" eff 



(11) 



From this we can convert our constraints on the turnover 
scale into constraints on fini/i 2 , N e g and subsequently 
z cq as follows. In Figure [4] we plot the ratio fim/fir as 
a function of the turnover scale (fco) measured from our 
CAMB power spectra for a range of N e ff between and 
10. Each line represents results for a wide range of values 
for fiM ft 2 , which increases with increasing fco- Projecting 
an observed turnover scale onto this sequence of curves 
allows us to map any value of fco to a unique set of values 
for fiM^i 2 and fiM/fir = l+z cq as a function of N e g. 

To see the form of the resulting constr aint in the 



TVcff — fiM ft 2 plane we can rearrange Equation 11 to get 



iVe ff = 4.403 



fiM 

fi7 



(fi M ^ 2 



(12) 



where the first bracketed factor is constrained from the 
turnover measurement and the second bracketed factor 
is known precisely from the CMB. Hence, we expect the 
scale of the turnover to place a roughly linear degenerate 
constraint in this plane, since our one measurement can 
not place closed constraints on two values. 

We also show on Figure [5] the WMAP CMB con- 
straint in the N e ff— flyih 2 plane. We made use of the 
ACDM+iVcfr WMAP 7-year chains, downloaded from the 
WMAP LAMBDA Data products siteQ and analysed 

t Note that due to the precision with which the CMB 
temp erature is measured (Tcmb=2.72548±0.00057K, |Fixsen| 
12009b , the photon energy density is precisely known to be 
n 7 /^=2.47274xl0- 5 . 

i http:/ /lambda. gsfc.nasa.gov/product/map/current/ 
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Figure 5. The constraints in the N^ff-Vlyih plane derived 
from the WiggleZ measurement of the turnover scale (in green) 
compared to CMB CosmoMC constraints from WMAP (in yel- 
low) . Dark shaded regions depict 68% confidence intervals and 
light shaded regions depict 95% confidence intervals. Grey con- 
tours depict lines of constant Qyi/Q T =l+z e q. 



them using the GetDist package that comes as part of 
CosmoMC ( |Lewis fc Bridle|20rjf| ) . 

Immediately obvious from this plot is the degener- 
acy in the WiggleZ constraints expressed by Equation 
|12| The WMAP constraints also exhibit a strong degen- 
eracy, although oriented differently and along the line 



l+^ C q~3200. While this degeneracy does not close for the 
WiggleZ constraint, it does for the WMAP constraint, 
demanding N e «>2. While we do not perform a rigorous 
joint fit of these measurements, it is clear that the Wig- 
gleZ turnover constraint should help to close the CMB 
constraints at an upper limit N c s<9. 

While our measurement of the WiggleZ turnover 
scale clearly provides little constraint on N e g (even with 
a strong prior on Q.Mh 2 ), a constraint on Qy[h 2 can 
be derived given a prior on N c a. Assuming the value 
N e g=3, we obtain a WiggleZ turnover constraint of 
n M h 2 =0.136t° ° f 2 (68%); ±£2™ (95%). As illustrated by 
Equation [TT] given a strong prior on AW, the fractional 
constraint on flMh 2 maps almost identically to a frac- 
tional constraint on z cq . Thus, we find that the WiggleZ 



turnover provides the constraint 



=3274 



+631 
-1259 



(68%); 



-i79i (95%). All results are unchanged if the expected 
value A^efi = 3.046 is taken for the assumed prior instead. 



4 FORECASTS FOR FUTURE SURVEYS 

From these results, the question naturally arises: what 
sort of constraint on N c ff, Q,mh 2 and z eq would be possi- 
ble from a larger survey such as BOSS or Euclid, encom- 
passing a larger volume and hence, providing a stronger 
turnover constraint? 

In this section we examine the prospects of using the 
turnover position in current and future surveys for cos- 
mology constraints. The key questions we seek to answer 
are: how effective will turnover constraints on N e s and 
QmIi 2 be from BOSS and Euclid and how big does a sur- 
vey need to be for these constraints to be competitive 
with those from CMB observations? 

We will find that the required volumes are large but it 
is important to note that photometric redshift surveys for 




Volume[(Gpc//i) 3 ] 



Figure 6. The precision of turnover scale (fco) measurements 
and subsequent constraints on (Jm/i 2 and JV e ff as a function of 
survey volume. Vertical dotted lines indicate the approximate 
volume of WiggleZ, BOSS and Euclid in order of increasing 
volume. For the SIm'' 2 constraints an A c g=3 prior is assumed 
and for the 7V e ff constraints an Qi^h 2 = 0.135 prior is assumed 
(although results are insensitive to these choices). Dark shaded 
regions depict 68% confidence intervals and light shaded re- 
gions depict 95% confidence intervals. 



which imaging systematics are under good control are just 
as effective for this science as spectroscopic redshift sur- 
veys. This is because the distance uncertainties associated 
with photometric redshift errors will be small compared 
to the scales of the turnover ( Blake fc Bridle|2005 1 . 



4.1 Constructing mock survey constraints 

We have repeated the analysis presented in Section [3] on 
a series of mock power spectra generated using CAMB 
in our fiducial cosmology and given ranges in k and noise 
properties designed to represent future surveys. We ignore 
the complicating issues of survey selection and geometry 
and assume that the scales relevant to this analysis remain 
volume limited in their precision. Under these conditions, 
survey volume (denoted as V) is the only survey param- 
eter relevant to our calculation. In all cases, we assume 
that the largest scale probed is L max = \fV and choose a 
power spectrum binning given by Afc = 27r/L max . In each 
case we combine 20 chains for the analysis of each fore- 
casted survey volume, each with differing random seeds 
for the noise generation process. 

For these calculations we have changed the maximum 
value of k over which we perform our turnover fits to 
fcmax=0.04 [ft/Mpc]. This was done to avoid a system- 
atic bias introduced by the BAO features of the power 
spectrum which becomes significant when volumes exceed 
those of WiggleZ. This bias is driven by a degeneracy be- 
tween /3 and feo and a dependance of j3 on fc max . 

We also slightly change the functional form of our 
model power spectrum for this calculation. With suf- 
ficiently large volume, the primordial power spectrum 
should begin to emerge in a survey's largest modes. The 
primordial power spectrum is expected to be a power law 
(or nearly-so), which differs substantially from the sim- 
ple asymmetric logarithmic parabola model of Equation 
[4] To capture this behaviour in our forecasts for surveys 
with very large volumes, we instead use for our power 
spectrum model an asymmetric logarithmic hyperbola on 
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0.12 0.14 0.16 0.18 



Figure 7. The constraints in the N e ff-£lMh 2 plane derived 
from our forecasts for the BOSS (in orange) and Euclid (in 
blue) measurements of the turnover scale. These are compared 
to the CMB constraint from WMAP (in yellow) and our Planck 
forecast (in grey). Dark shaded regions depict 68% confidence 
intervals and light shaded regions depict 95% confidence in- 
tervals. Grey contours depict lines of constant Q]^i/Q T =l+z eq . 



scales larger than the turnover. Equation[4]then becomes: 
log 10 P +p-^/a 2 x 2 +p 2 if k < ko, 



log 10 P(fe) 



log 10 P (l- fix 2 



if k ^ ko- 



(13) 

This model adds an additional degree of freedom (p) 
which describes how quickly P(k) asymptotes to a power 
law at k<kQ. We have verified that using this model on 
the observed WiggleZ P(k) and on model power spectra 
of comparable volume generates consistent results for the 
turnover scale. 



4.2 Forecast results 

The results of our survey forecasts are shown in Figure 
[6] where we present, as a function of survey volume, the 
precision of the turnover scale measurement and of the 
resulting constraints on Qyih 2 and A^ff. For the Om/i 2 
constraints an N c b=3 prior is assumed and for the iVoff 
constraints an Qyih 2 =0.135 prior is assumed. Results are 
insensitive to these choices however, particularly for large 
volumes. 

In this figure we highlight the results for survey vol- 
umes approximately equal to those of WiggleZ, BOSS 
and Euclid representing a series covering the range of 
our calculations in logarithmically equal steps in volume. 
The BOSS survey ( [Eisenstein et al.|2011[ ) will map 10000 
square degrees, collecting spectroscopic redshifts for lumi- 
nous red galaxies out to z~0.7. With completion planned 
for 2014, BOSS represents the largest survey for which 
near-term results will be possible. The Euclid satellite 



(Laureijs et al. 20111 will conduct a wide-field (15000 



square degree) extragalactic survey covering redshifts out 
to z~2. With a target launch date of 2017-2018, the an- 
ticipated specifications of this survey represent a rough 
limit to our capabilities for conducting the science dis- 
cussed here at the turn of the next decade. 

In Figure [fi] we see that results for V~l [(Gpc/h) 3 ] 
are roughly consistent with the WiggleZ constraints (al- 
though slightly better, presumably due to the lack of co- 
variance in our forecast model) presented in Figures H] 



and [5] a k precision of ±20% (68%); ±40% (95%) with 
little resulting constraint on iV c ff<10. The 24% precision 
of the constraint on Q,yih 2 (and thus, z cq ) for this case 
maps almost directly to the precision of the constraint on 
ko- 

As volume approaches that of BOSS at V~10 
[(Gpc/h) 3 ], the nature of these constraints changes signif- 
icantly. With a sizeable increase in the number and accu- 
racy of measured modes and scales beyond the turnover, 
the constraints on ko improve dramatically to a precision 
of ±9% (68%); ±16% (95%). This translates to signifi- 
cantly improved constraints on Q,yih 2 and z eq of ±10% 
(68%); ±20% (95%). The constraints on N eS start to be- 
come interesting at a level of ±5*% (68%); ±gf a % (95%), 
as well. 

Looking past BOSS to much larger future surveys, 
precision improves more slowly with increasing volume. 
However .with V / ~100 [(Gpc/h) 3 ], the proposed volume 
of EuclioTj is a dramatic order-of-magnitude increase over 
BOSS. As a result, Euclid should ultimately be able to 
measure k to a precision of ±3% (68%); ±5% (95%) 
with corresponding constraints on Vlyih 2 and z cq of ±4% 
(68%); ±6% (95%) and N cS of ±%% (68%); ±f 5 % (95%). 

We contrast these forecasts with the results antici- 
pated from Planck. For our Planck constraints we gener- 
ated simulated Planck data following the procedure used 
111 |Perotto et al] |2006l ) using the specifications for the HFI 
bolometers given in the Planck Blue Book ( Collabora- 
tion" |2006 1. We then used CosmoMC to generate chains, 
assuming a ACDM+A'ofj model. 

In Figure [7] we illustrate our forecasted constraints 
for BOSS and Euclid in the N e g-QMh 2 plane, compared 
to the CMB constraints of WMAP and our forecasts for 
the Planck Surveyor. We can see from this figure that 
BOSS should make important contributions to the present 
WMAP constraint, closing the upper limit of the 68% 
confidence contours at AT e g<6. Looking to the future how- 
ever, we see the dramatic improvement Planck will make 
to this measurement. In detail, we find the constraints 
from Planck on (S7m/i 2 , N e g, z cq ) should be ±3, ±9, and 
±2% (68% confidence); ±5, ±14, and ±4% (95% confi- 
dence) respectively. While this constraint on z cq is signif- 
icantly better than what Euclid will provide on its own 
(even with strong priors), the constraints on Qm/i 2 and 
A^eff are otherwise comparable. 



5 SUMMARY AND CONCLUSIONS 

We have presented an analysis of the WiggleZ Dark En- 
ergy Survey, constructing a galaxy power spectrum op- 
timised for studying the largest scales of the survey. We 
extract from this the most robust measurement to date of 
the scale of the turnover in the Universe's matter power 
spectrum. From this, we have obtained the first distance 
measurement and cosmology constraints yet derived from 
measurements of this feature. We have also constructed 
forecasts for the precision of this analysis for future sur- 
veys, contrasting the constraints we expect to obtain from 
turnover measurements in future redshift surveys to com- 
plementary constraints from the published observations 
of WMAP and our forecasts for the Planck Surveyor. 
Details of our results from analysis of the WiggleZ 



" From the Euclid Science Requirements Document, ESA Sci- 
ence document reference number DEM-SA-Dc-00001. 
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Table 2. A summary of our WiggleZ turnover-derived cosmology measurements and of our forecasts for the precisions of turnover- 
derived and CMB results in future/ongoing surveys. Unbrackctcd quantities are 68% confidence results and bracketed quantities 
are 95% confidence results. 



Survey 


ko 








2cq 








WiggleZ Measurement 


o nifin+ 0035 

U.U1DU_q 0041 


(-+0.0073-, 
*.-0.0075-> 


U.1JD_ 052 


f+0. 073-, 
V-0.074J 


397A+ 632 
JZ '*-1260 


,- + 1757-, 
( .-1791^ 






BOSS Forecast Precision 
Euclid Forecast Precision 
Planck Forecast Precision 


±9% 
±3% 


(±16%) 

(±5%) 


±10% 
±4% 
±3% 


(±20%) 

(±6%) 
(±5%) 


±10% 
±4% 
±2% 


(±20%) 
(±6%) 
(±4%) 


-56 /o 
-21 /o 
±9% 


( + -\f%) 

(±35%) 

(±14%) 



dataset (all uncertainties quoted at 68% confidence) are 
as follows: 

• We present an MCMC method which removes win- 
dow function convolution effects while coadding the indi- 
vidual observed power spectra of a galaxy survey observed 
over several disconnected regions. 

• Applying this method to WiggleZ data we find that 
the survey is able to probe modes at-and-beyond the scale 
of the turnover. 

• Using an MCMC approach, we have fit an asym- 
metric logarithmic parabola to the observed WiggleZ 
power spectrum. From this analysis we find the scale of 
the turnover to be k =0.0160to%oli- Tnis is in excel- 
lent agreement with our fiducial standard ACDM value 
of k = 0.016 [Mpc/ft]. 

• Parameterising the power spectrum beyond the scale 
of the turnover as P(k<ko)ock n , we find n> — 1 at nearly 
95% confidence. The standard ACDM value of n=l is 
easily accommodated by our fit. 

• The continuance at large scales of the small-scale 
asymptotic value of n—— 3 is completely ruled-out by our 
analysis. 

• We have performed the first measurement of the peak 
position in the cosmological power spectrum, representing 
the first secure and quantified observation of this feature 
to date. 

From this measurement we then extract - for the 
first time using a measurement of the turnover scale - the 
following information: 

• A model-independent distance measurement to 
2=0.62 in units of the z=3145 horizon scale at 
the redshift of matter-radiation equality (m) of 
Dv(zeff=0.62)/r H =18.3t|i 

• A measurement of the cosmological density parame- 
ter r2M/i 2 =0.136l U Qg2 (assuming a JV ff=3 prior). 

• A measurement of the redshift of matter-radiation 
equality z cq =3274ti26o (assuming a N c g—S prior). 

Looking to the future, we have computed forecasts 
for the turnover precision attainable by BOSS and Eu- 
clid. We find that BOSS should substantially improve 
upon the results presented here, reaching precisions in 



standard ACDM predictions for P(k) on scales compara- 



ble to and larger than the turnover (Dalai et al 



2008 



(ko, Q,M.h 



,N eS ) of (±9, ±10, ±10, 



This repre- 



sents sufficient precision to sharpen the constraints on 
N c fj from WMAP, particularly in its upper limit. For 
Euclid, we find corresponding attainable precisions of 
(±3, ±4, ±4, ^L^i)%- This represents a precision approach- 
ing our forecasts for Planck. 

We emphasise that these results are all obtained 
within the theoretical framework of Gaussian primor- 
dial fluctuations. In the event that cosmological struc- 
ture formation was seeded by non-Gaussian fluctuations, 
scale-dependent bias effects may substantially change the 



Dcsjacqu es et aL"1|2009[ |2011b|a| |Giannantonio fc Por-| 



ciani|2010 1 . Measurements of primordial Gaussianity from 



studies of large-scale structure have proven to provide 
constraints of comparable precision to those from the 
CMB ( Slosar et al.|2008" I with forecasts for future surveys 
suggesting that the two approaches will remain comple- 
mentary for some time ( Giannantonio et al. 20121. In- 



terestingly, current constraints allow for cases where the 
power spectrum exhibits no turnover but merely an inflec- 
tion at turnover scales, with P(k) being a declining func- 
tion of k at large scales. Many such models are compati- 
ble with the results presented in this study but the lower 
limit we place on n in the large-scale limit of P(k)ock n 
suggests that the WiggleZ power spectrum measurement 
may be capable of constraining their allowed parameter 
ranges. However, such effects are expected to scale with 
galaxy bias b as (6—1) {e.g. Dalai et al.||2008 l and Wig- 
gleZ galaxies are known to have a bias near unity (e.g. 
Bla ke et al.||2010| ), greatly reducing the degree of such 
effects for this survey. Regardless, we forgo this analysis 
for now and focus instead on the case of purely Gaussian 
primordial fluctuations. Future studies involving WiggleZ 
bispectrum measurements will address this science. 

Furthermore, while this paper focuses on the scale of 
the turnover as a cosmological constraint, there is much 
more information present in the full shape of the power 
spectrum which could provide (for example) constraints 
on JIm/i 2 (see Parkinson et ah, submitted), tightening the 
constraints presented here not only on this parameter but 
on N e g as well. Such a situation is similar to that of BAO 
studies where final cosmological constraints tend to in- 
corporate information from both the BAO scale and the 
shape of the underlying power spectrum (see Figure 6 of 
Blake et al.|20lfb for example). As such, the constraints 
presented here likely represent underestimates of the true 
potential of these experiments. 

Lastly, with increased constraints on N e g and Q,yih 2 , 
the absolute scale of the turnover could be accurately cali- 
brated, permitting its use as a standard ruler for measur- 
ing the distance-redshift relation. Our forecasts suggest 
that a regular volume of V~10 [(Gpc/h) 3 ] with a well 
defined effective redshift would permit a distance mea- 
surement of roughly 10% accuracy. Euclid for instance 
could provide several such volumes arranged across red- 
shift, perhaps enabling the first measurement of a Hub- 
ble diagram using the turnover scale as a standard ruler. 
Given the purely linear evolution of the matter power 
spectrum on scales of the turnover, this could provide a 
powerful check against systematic redshift-dependent bi- 
ases in BAO studies. 
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